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Abstract 

We report a variational approach to the nonlinearly screened interaction of charged particles 
with a many-electron system. This approach has been developed by introducing a modification 
of the Schwinger variational principle of scattering theory, which allows to obtain nonperturbative 
scattering cross-sections of moving projectiles from the knowledge of the linear and quadratic 
density-response functions of the target. Our theory is illustrated with a calculation of the energy 
loss per unit path length of slow antiprotons moving in a uniform electron gas, which shows good 
agreement with a fully nonlinear self-consistent Hartree calculation. Since available self-consistent 
calculations are restricted to low heavy-projectile velocities, we expect our theory to have novel 
applications to a variety of processes where nonlinear screening plays an important role. 
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The interaction of external charges with solid targets is a problem of both fundamental 
and practical interest in contemporary physics.—^ Since the targets of interest, either bulk 
solids, surfaces, or structures of lower dimension typically consist of many interacting elec- 
trons, a remarkable progress has been achieved along the lines of the so-called 'dielectric' 
perturbative formulation of scattering. This formulation utilizes density-response functions 
of many-body systems to characterize their dynamical reaction to external perturbations.^ 

Linear-response theory provides a qualitative description of the energy loss per unit path 
length of external particles of charge Z\ interacting with solid targets, i.e., the so-called 
stopping power (SP) of the solid. However, it yields SP that is proportional to Z\ and 
cannot, therefore, account for the existing differences between the ranges of positive and 
negative pions^ and the slowing of protons and antiprotons.^ The quadratic-response treat- 
ment considerably improves the description of the interaction of external charges with solid 
targets,^ and yields a Z\ correction to the energy loss that accurately accounts for the 
measured energy loss of protons and antiprotons in the high-velocity regime?*^ Consider- 
able progress has been achieved recently in the framework of quadratic-response theory; n > 12 
however, at low velocities this theory is only able to quantitatively account for the strong 
influence of unit-charge projectiles {Z\ = ±1) in the case of high-density targets. 

Another approach to investigate the interaction of external charges with a many-electron 
system is based on the ordinary formulation of potential-scattering theory. In this ap- 
proach , 13 i 14 i 15 i 16 i 17 SP for a heavy particle is determined in the low-velocity limit from the 
knowledge of the scattering phase shifts corresponding to a static nonlinearly screened po- 
tential. These potential-scattering calculations include all orders in the projectile charge; 
however, they have the limitation of being restricted to low velocities {y « vp, vp being 
the Fermi velocity) of heavy projectiles moving in bulk materials. 

In the case of few-body scattering, variational methods of the Schwinger type are known 
to provide nonperturbative representations of the collision matrix, 18 which in the case of the 
scattering of external charges would not be restricted to low velocities and heavy projectiles. 
However, the direct application of these methods to solid targets composed with a large 
number of electrons is not feasible, since it requires the knowledge of the many-particle 
ground and excited states of the target. 

In this paper, we put forward a theory which reconciles the dielectric formulation of scat- 
tering with the variational method. To do this, we introduce a modification of the Schwinger 
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variational principle of scattering theory to provide a nonperturbative representation of the 
scattering cross-section of external charges interacting with a many-electron system. The 
novelty of this method is that it does not require the knowledge of the many-particle ground 
and excited states of the target; instead, only the linear and quadratic density-response func- 
tions are needed, which have been previously obtained in the framework of time- dependent 
density-functional theory (TDDFT).— 

Because of the quite general nature of this approach, it can be readily applied to the 
investigation of a variety of processes involving the inelastic scattering of charged particles 
by many-electron targets, such as the SP for moving ions, electron and positron energy-loss 
spectroscopy, inelastic low-energy electron diffraction, and hot-electron dynamics. As our 
formula only requires the knowledge of the two lowest-order density-response functions of 
the target, the implementation of this approach has a computational cost equal to that of 
available quadratic-response theories. 

Let us consider the scattering of a charge by an arbitrary many-electron target, in which 
the target is scattered between its ground |0) and excited \n) states. We assume that 
the initial and final states = |p i; 0) and tyf = \pf,n) of the projectile-target interacting 
system contain a free particle of charge Zi and momentum pj and pj, respectively. According 
to the bilinear form of the Schwinger variational principle,— the functional 

[T fi ] = (^\V\^,) + ^ f \V\^i)-(^\V~VGi +) V\^t) (1) 

is stationary under the variation of the trial (unknown) eigenstates \I/jr and ^ of the full 
interacting hamiltonian and at its stationary point gives the exact transition-matrix elements 
between the (known) initial and final tyf states. Here V is the Coulomb interaction 
between the projectile and the target, and is the Green function associated to the 
hamiltonian of the projectile-target system without the mutual interaction V. 

The fractional form of the Schwinger variational principle is obtained by substituting (see 
Ref. Q, p. 412) 

¥+ - A¥+ *7 - B^J (2) 

in Eq. and treating the coefficients A and B as variational parameters. If one then 
returns to Eq. ([!} and uses the free states for the trial functions, one finds 

Tfi = Tff/ (l - if /if) , (3) 
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where TL' = (^f\V\^i) and Tj { = f\VG V|^) represent the first two Born transition 
amplitudes. 

Equation (J3J) has proven very useful in atomic scattering, where the |0) and \n) states 
of the target can be known, at least approximately. However, in the case of solid targets 
the many-body ground and excited states are difficult to know. In order to find a represen- 
tation of the scattering cross-section of external charges interacting with a many-electron 
system in terms of density-response functions, which can be calculated from the knowledge 
of one-electron states, we first construct the differential cross-section for the projectile to 
be scattered between the free-particle states of momenta p and p — k (we use atomic units 
throughout) 18 

da/dk= (16n 4 /v) ^2\T fri:i0 \ 2 5(uj - u n0 ), (4) 

n 

where the sum is extended over a complete set of eigenstates of the target, u n0 are the 
excitation energies of the target, uo = k- v — k 2 /2M, M is the projectile mass, and v = p/M 
is the projectile initial velocity. We then introduce the functional 

[da/dk] = (16tt» O*/^**)) + (*fn\V\*to) 

n 

-(^ n \V - VG^V\*+)\ 2 5{u - oo n0 ), (5) 

we compare this functional with that of Eq. ([Tfl. and noting that a linear combination of the 
squares of the absolute value of a stationary quantity is a stationary quantity, we conclude 
that the functional of Eq. (JSJ) gives the exact cross-section of Eq. at its stationary point 
with respect to the variation of the trial states ^f and ^J n - The variational principle based 
upon the functional of Eq. (|SJ) is the 'global' analog (summed over all final states of the 
target) of the corresponding 'local' variational principle based on the functional of Eq. 
for the transition to a particular state. These two principles are both exact and equivalent 
to each other. 

In order to obtain the fractional form of the global variational principle, we first apply 
the substitution (J2J) to the functional and we then approximate the trial functions by 
the free states. Hence, the problem reduces to finding the stationary value of the functional 

[da/dk] = 16ir 4 /v x 

£ | (A + B - A B)T^ + AB 5(u - u n0 ) (6) 

n 
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with respect to the variational parameters A and B. Equation (J6J) can be rewritten as 

[da/dk] = (A + B- AB) 2 R 2 + 

AB(A + B- AB)R 3 + A 2 B 2 R 4 . (7) 

Here, R 2 , R3, and R4 represent contributions to the differential cross-section that are pro- 
portional to Z 2 , Zf, and Zf, respectively, which can be obtained from the knowledge of the 
corresponding density-response functions of the target, as shown in Ref. Il2. 

In what follows, we neglect the last term of Eq. (J2J) (see the discussion below). Solving 
the variational problem for the function of Eq. ((7J) with respect to the parameters A and B, 
we find its stationary value as 

da/dk = R 2 f(R 3 /R 2 ), (8) 

where 



/(*) 



16 - 32a; - 56a; 2 + 72x d - 27x 4 



+ (2 - x) (4 - Ax + 9a; 2 ) 3/2 ] / 32 (1 - xf. (9) 

Equation (|SJ) is the 'global' analog of the 'local' Eq. (jSJ). If the quantity x = R 3 /R 2 is small, 
we may expand the expression Q and obtain 

da/dk = R 2 + R 3 + . . . , (10) 

which agrees with the Born series through second order in the interaction between the 
projectile and the target. 

In the case of antiprotons {Z\ = — 1) moving slowly in a metal, the Zf contribution 
to the scattering cross-section R 3 , which is negative, can overcome for low-density metallic 
targets the positive Z\ contribution R 2 , thereby leading to a physically meaningless negative 
scattering cross-section. In contrast, our nonperturbative differential cross-section (jHJ) is 
positively defined, as shown in Fig. [T] where the function f(x) of Eq. Q is plotted. 

We illustrate our theory with a calculation of the SP [—(dE/dx)} for slow protons and 
antiprotons moving in a uniform electron gas (EG) of density n characterized by the density 
parameter r s = (3/47r)n . At low velocities the projectile-target Coulomb interaction 
is relatively strong, so this represents an unfavorable situation for a linear or quadratic 
perturbative approach. 



FIG. 1: f(x) of Eq. © (solid line). If x is small, f(x) = 1 + x + . . . (dashed line), which yields 

Eq. dini). 



The SP is obtained by multiplying the differential cross-section (da/dk) by the energy 
transfer u and integrating over the momentum transfer k. As the mass M of our projectile 
is much larger than the electron mass, uo ~ k • v and one writes 



where Vt is the normalization volume. We have evaluated the differential cross-section of 
Eq. (jSJ) from the knowledge of the linear and quadratic density-response functions of the 
uniform EG, which we have calculated in the random-phase approximation (RPA).— For 
slow projectiles {y < vp), the SP is found to be proportional to the projectile velocity v. 

In Fig. |21 the variational SP of a uniform EG for slow protons and antiprotons (thick solid 
lines) is compared to the corresponding quadratic SP [obtained from Eq. ()10|)] (dashed lines) 
and accurate fully nonlinear potential-scattering calculations that we have carried out along 
the lines of Refs. Q and but with exchange-correlation (xc) excluded (chained lines). 
In the high-density limit (r s — > 0), the Born series quickly converges and all calculations 
coincide. However, as the electron density decreases, a correct description of nonlinear 
interactions requires to go beyond quadratic-response theory. 

In the case of antiprotons, quadratic-response calculations (lower dashed line of Fig. |2J) 
overestimate the negative nonlinear contribution to the SP (see also Fig. Q), which for 
r s > 5.5 becomes larger in magnitude than the linear (Zf) term. In contrast, the variational 
SP (lower thick solid line) is positive for all electron densities and shows good agreement 
with the fully nonlinear potential-scattering calculation that we have performed by solving 
self-consistently the Hartree equation of a static antiproton (chained line with circles). Since 
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the Hartree SP is obtained from accurate phase-shift calculations to all orders in Zi, the 
agreement between variational and Hartree calculations gives us confidence in the accurate- 
ness of the variational approach. 

In the case of protons, our variational SP (upper thick solid line of Fig. |2J overestimates 
the potential-scattering self-consistent Hartree calculations (chained line with squares) even 
more than within quadratic-response theory (upper dashed line). This is due to the fact 
that in its present (initial) form our theory is only applicable to point charges which do not 
support bound states or resonances.— Protons in a uniform EG are known to support either 
a bound state (at r s > 2) or a resonance below the Fermi level,— which can be inferred by 
the behavior of our calculated scattering phase shifts. Electrons accommodated in either 
a bound state or a resonance efficiently screen the projectile, thereby decreasing the SP as 
shown in Fig. EJ 

The first unambiguous evidence for a velocity-proportional electronic SP of solid targets 
has been reported recently, with measurements of the energy loss of slow antiprotons in Ni, 
Au, C, and Al.^ At low velocities, the energy loss of ions in metals is mainly due to the 
slowing by valence electrons, which can be approximately described with the use of a uniform 
EG. 20 Hence, we have also plotted in Fig. ED the friction coefficients —{dE/dx)/v obtained 
from the slow-antiproton energy-loss measurements of Refs. |5j and [g, at r s =1.20, 1.35, 
1.53, and 2.07 corresponding to the average density of valence 3<i 8 4s 2 , 5d 10 6s 1 , 2s 2 2p 2 , and 
3s 2 3p 1 electrons in Ni, Au, C, and Al, respectively. The variational and Hartree calculations 
for antiprotons are considerably closer to experiment (solid circles) than their quadratic 
counterpart. Our variational calculations for protons are also close to experiment; however, 
this might be due to a cancelation of the impact of bound and resonance states with the effect 



13, 3, and 
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of xc not included in our approach. The xc effects were included in Refs. 
framework of density-functional theory (DFT),— and they can be included in our variational 
procedure by going beyond the RPA in the description of the density-response functions of 
the solid. Work in this direction utilizing the recently obtained quadratic density-response 
function from TDDFT with xc included 22 is now in progress. 

Although an additional approximation of omitting the last term in Eq. (JJJ) has been made, 
we note that by including it a more general analytical formula for the differential cross- 
section would be obtained. However, while the terms R<± and R3 require the knowledge of 
the linear and quadratic density-response functions, respectively, which are available from 
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FIG. 2: Variational (solid lines) and quadratic (dashed lines) SP of a uniform EG for slow antipro- 
tons (lower curves) and protons (upper curves), divided by the projectile velocity, as a function 
of r s . The chained lines with circles (squares) represent fully nonlinear Hartree calculations for 
antiprotons (protons), performed along the lines of Refs. [lj and Il5l but with xc excluded. The 
circles (squares) show the experimental results for antiprotons (protons) reported in Refs. QandO. 



the literature for a number of systems of physical importance, R4 is related to the cubic 
response function, which has not been determined as yet. Although the generalization of 
the present theory in this direction is conceptually straightforward, for the sake of simplicity 
(and implementability in calculations at present) we defer it until later publications. 

In conclusion, we have reported a new nonperturbative variational approach to the non- 
linearly screened interaction of charged particles with a many-electron system, which goes 
beyond the conventional linear and quadratic theories. This approach has been developed by 
introducing a modification of the Schwinger variational principle of scattering theory, which 
allows to obtain nonperturbative scattering cross-sections from the knowledge of the linear 
and quadratic density- response functions of the target. Our approach, which includes contri- 
butions to all orders in the projectile-target Coulomb interaction and agrees with the Born 
series through second order, represents a considerable improvement over the quadratic the- 
ory (Zf) approximation; in particular, our variational differential cross-section is positively 
denned, which is known not to be the case for the quadratic theory. 

We have illustrated our theory with a calculation of the stopping power of uniform EG 
for slow protons and antiprotons, which in the case of antiprotons shows good agreement 
with fully nonlinear Hartree calculations. Our calculations indicate that by going beyond 
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quadratic theory the variational procedure considerably improves the agreement with recent 
measurements of the stopping power of Ni, C, Al, and Au for slow antiprotons, though xc 
effects still need to be taken into account. In the case of protons, the presence of bound 
states and resonances also needs to be incorporated into our theory. Since self-consistent 
Hartree calculations have the limitation of being restricted to low heavy-projectile velocities 
(v << vf), we expect our theory to have novel applications in the investigation of a variety 
of processes involving the inelastic scattering of charges by many-electron targets. 
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